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Abstract 

A family of multispecies drop-push system on a one-dimensional lattice is 
investigated. It is shown that this family is solvable in the sense of the 
Bethe ansatz, provided a nonspectral matrix equation is satisfied. The 
large-time behavior of the conditional probabilities, and the dynamics of 
the particle-type change are also investigated. 
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1 Introduction 



Various aspects of one-dimensional asymmetric exclusion processes have been of 
physical interest. These contain, for example, the kinetics of biopolimerization 

dynamical models of interface growth 2,^ ^^^d the traffic models This 
model is also related to the noisy Burgers equation 0] , and hence to the study of 
shocks I^IIHI ■ The dynamical properties of this model have also been extensively 
studied, for example in [H1[71IH1- 

In the study of stochastic processes, the term solvable has been used in 
several meanings. In [5l ll(Jl[mil2[ll3l[T^ll5M 16' . solvability means solvability 
in the sense of the Bethe ansatz, or factorization of the A^-particle scattering 
matrix to the two-particle scattering matrices. This is related to the fact that for 
systems solvable in this sense, there are a large number of conserved quantities. 
In | T7^18.19-.20..21 22 23 24 25 ■ 26; , solvability means closedness of the evolution 
equation of the empty intervals (or their generalization). And in |27ll28ll^ . 
solvability means that the evolution equation for the n-point functions, contain 
only n- or less- point functions. 

In 1^, the Bethe ansatz was used to solve the asymmetric simple exclusion 
process on a one-dimensional lattice. In |10| . a similar technique was used to 
solve the drop-push model ,3Q., and a one-parameter family of reactions con- 
taining the simple exclusion processes and the drop-push model as special cases. 
In |12| . the same technique was used to solve a two-parameter family of pro- 
cesses, involving bidirectional diffusion, exclusion, and pushing. The behavior 
of this last model on the continuum, was investigated in In all of the above 
cases, the essence of the method has been to replace the reactions by suitable 
boundary conditions. The boundary conditions involved were extended in |15| 
to describe systems in them the process of annihilation exists as well, and in 
to general boundary conditions in the continuum. 

All of the above studies have been about single-species systems. In |13II14) . 
systems with exclusion processes were investigated, which contained more than 
one species. It was shown in J^], that in order that such systems be solvable 
in the sense of the Bethe ansatz, certain relations should be satisfied between 
the rates. This relations can be written as some kind of a spectral Yang-Baxter 
equation. In |14j . it was shown that this spectral equation is equivalent to a 
nonspectral matrix equation involving the rates. 

Here we want to extend this approach to the case of drop-push models. The 
systems under consideration, consist of N particles, which can be of several 
species. This particles live on an infinite one-dimensional lattice, so that each 
site of the lattice is either empty or contains one particle. Each particle hops to 
the site at its right-hand side with unit rate, if that site is empty. If that site is 
occupied, then the particle may push the other particle, and at the same time a 
reaction may occur between these neighboring particles changing their species. 

The scheme of the paper is the following. In section 2, a multi-species exten- 
sion of the drop-push model is introduced and the use of a suitable boundary 
condition instead of the reaction is investigated. In section 3, the solvability 
condition (in the sense of the Bethe-ansatz) for this reaction is investigated, 
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and it is shown that this condition, which is a spectral equation for the matrix 
of the reaction rates, can be rewritten as a nonspectral equation for the same 
matrix. In section 4, the conditional probability is obtained, and its behavior 
in the two-particle sector, specially its large-time behavior, is investigated. In 
section 5, the dynamics of the particle- type change in the two-particle sector 
is investigated and the large-time limit of the probability of particle-types is 
obtained. Finally, section 6 is devoted to the concluding remarks. 



2 Multi-species extension of the drop-push model 

In the ordinary drop-push model, the system consists of a single type of particles, 
living on a one-dimensional lattice. Each site of the lattice is either empty or 
occupied by one particle. Any particle can hop to the site at its right neighbor, 
with the rate 1, if that site is empty. If the right neighbor site is occupied, the 
particle can still hop to that site and push the second particle, with the same 
rate 1. One can write the reactions like 

AA- ■ ■ A9 (/}AA- ■ ■ A, with the rate 1. (1) 



In this reaction, the total number of particles is conserved. For a system con- 
taining N particles, the question of interest is to determine the probability of 
finding the N particles in sites xi to xn, where 

x.i <Xj, for i < j, (2) 

the so called physical region. It is easily seen that the evolution equation for 
this probability, . . . , x^; t), is 

Pixi, . . . ,XN\t) =P{xi - 1, . . ■,XN;t) H h P(xi, . . . ,XN -l;t) 

-NP{xi,...,XN\t), (3) 

if among the sites Xi, no two are adjacent; that is, if Xi < Xi^i — 1. For a block 
of {n + 1) adjacent sites, the evolution equation becomes 

P{xq — X, . . . ,Xn = X + n;t) =P{xq ~ x — 1 , . . . , Xn ^ x + n; t) + ■ ■ ■ 

+ P{xo — X — 1, . . . ,Xk = X + k — 1, 

Xk+i = X + k + 1, . . . ,Xn — X + n;t) 
+ {n + l)Pixo...,Xn;t). (4) 

This looks like different from (|2Jl. However, defining a boundary condition 

P{...,x,x,...):=P{...,x-l,x,...), (5) 

makes the forms of and (O similar. One notes that P{. . . ,x,x, . . .) has in 
fact no physical meaning, since the argument of P in that expression is not in 
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the physical region. But its introduction helps to solve the evolution equation, 
as it was done in j^lj (and in |2| for the exclusion process). 

Now suppose that the system consists of k species of particles; that is, an 
occupied site may have k different states. Assume moreover, that if the right 
neighbor of a particle is free, the reaction is the same as ordinary drop-push 
model, without changing the type of the particle, but there is a difference when 
two particles are adjacent to each other: the left particle does push the right 
one with unit rate, but in the mean time there is a probability that the types 
of the particles change. So we have reactions like 

— > 0Aq, with the rate 1, 
9AjAs with the rate f'^a/::- (6) 

Consider a consisting of two particles, and denote the probability that the first 
particle be at the site x and of the type Aa and the second particle be at the 
site y and of the type Ap, by P°'f^{x,y). The evolution equations become 

y; t) ^P^P{x - 1, y- 1) + P"^(x, y - 1; t) - 2P'-P{x, y-t), x<y-l, 
P°'l^{x, x + l;t) =P"'3(x -l,x + l:t) + b^^^^s P'^^x ~l,x]t) 

- P'^'^ix, x + l;t)- P°'P{x,x + 1; t), (7) 

where 

B'^P ■.= Y,h-^'^0. (8) 

7(5 

In fact, B"'^ is the overall pushing rate, in which the type change is unimportant. 
If this overall pushing rate is 1, the second equation in (|7|) is simplified and it 
is seen that it can be rewritten in the form of the first equation, provided one 
introduces the boundary condition 

P-^f^ix, X] t) = b'^'^-^s P^'Hx - 1, x; t). (9) 

One notes that all of the elements of the matrix b (including the diagonal ele- 
ments) are nonnegative, as they are rates, and b satisfies 

(s (g) s)6 = s ® s, (10) 

where 

s„ 1. (11) 

(This simply means that the sum of the elements each of the columns of b is 
equal to one.) A similar matrix b was introduced in |14| . however the diagonal 
elements of that b were not necessarily nonnegative. 

Now consider a system consisting of N particles of various species, with the 
evolution equation 

P(a;i, ...,XN;t) =P(a;i - 1, . . . , xat; i) H h P(xi, . . . , xa? - 1; t) 

~ NP{xi,...,XN;t), (12) 
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in the whole physical region, and the boundary condition 

P(. . . , Xk = X, Xk+i =x,...) := bk,k+i P(- • ■ , Xk = X - 1, Xk+i =x,...), (13) 
where 

6fc,fc+i 1 ® • • • (g) 1 ® (8)1 (g) • • • (g) 1, (14) 

and P is an A^-tensor the components of which are probabilities. It is seen that 
in this system, apart from the simple diffusion, there is a reaction between a 
block of n + 1 adjacent particles: 

■ ■ ■ Aa^% 0A^„ • • • A^„, with the rate (6„_i,„ • • • 6o,i)''°'"^"ao---a„- 

(15) 

This comes from the fact that 

P(a;o = X, . . . , Xn-i = X + n — l,Xn = x + n — 1) 
=ibn-i,n ■ ■ ■ bo,i)F{xo = a; - 1, ... , x^-i = x + n - 2, x^ = x + n - 1). (16) 

Note the order of the matrices b. This order suggests that if a collection of n + 1 
particles are adjacent, there is a probability that the first particle pushes the 
second and changes the type of the second (and itself) and then it is the second 
(modified) particle that interacts with the third. 

3 Solvability and the Bethe-ansatz solution 

Consider the evolution equation (|12|l with the boundary condition H13|) . To 
solve this equation, one as usual seeks the eigenvectors of the operator acting 
at the right-hand side of 112() . that is, one tries to solve 

E^{xi,...,xn) =*(a;i - 1,.. .,a;jv) + • •• + *(xi, . . . , Xjv - 1) 

-iV*(xi,...,xjv), (17) 

with 

*(..., Xfc = x, Xk+i = X, . . . ) := bk,k+i *(. . . , Xfe = a: - 1, Xk+i = x, . . . ). (18) 
The Bethe-ansatz solution to this equation is 

*(f) = ^A,e*"(P)"H, (19) 

where H is an arbitrary vector and the summation runs over the elements of the 
permutation group of N objects. Putting this is H17(l . one arrives at 

N 

iJ = ^(e-^"^-l), (20) 

k=l 
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while ((TH|l gives 

[1 _ (f'^) bk,k+i]A„ + [1 - (^"=+1) bk,k+i]A„a, = 0, (21) 

where 

{Pk+i, 3 = k 
Pk, + 1 . (22) 

From H21|) . one arrives at 

^crcTfc = 5'fe.fe+l[(7(pfc)ia-(P/c+l)]^o-, (23) 

where 

Sk,i{Pi,P]) ■= -(1 - Zj &fej)"^(l - &fe,i): (24) 

and 

. (25) 



So one can construct A^^s from Ai, by writing tr as a product of fr^'s; that 
is, one can write the A^-particle scattering matrix A, as a product of the two- 
particle scattering matrix S. However, as the generators of the permutation 
group satisfy 

CTfeCTfe+lCTfe = O-fc+lffcCrfc+l, (26) 

one also needs 

AakCTk+lCk ~ ^o"fc + icrfe(Tfc+i • (27) 

This, in terms of the S becomes 

5'l,2(P2,P3) S2,3{Pl,P3) Sl^2{Pl,P2) = 5'2,3(Pl,P2) 5*1,2(^1, Ps) 5'2,3 (P2 , P3 ) • (28) 

In terms of the i?-matrix defined through 

Sk,k+1 =■ ^k,k+l Rk,k+1, (29) 

becomes 

^2.3(P2,P3) RiAPI^Ps) Rl.2{Pl,P2) = 2(^1,^2) RiAPI^Ps) ^2.3(^2,^3), 

(30) 

which is the spectral Yang-Baxter equation. 

The Bethe-ansatz solution exists, iff the scattering matrix satisfies (|28|l . This 
is a general statement not coming from a specific reaction. In the problem 
studied here, S is of the form (|24|l . Comparing this with the 5- matrix obtained 
for the multi-species simple exclusion process |14j . it is seen that S{pi,pj) in 
these two different problems are transformed to each other by a simple change 
Zi <-> Zj. The definition of z in terms of p in the present paper is, however, 
different from that of jj^. It is also seen that with the changes 61^2 ^ ^2,3 and 
zi ^ Z3, (|28|l is transformed to the spectral equation for S obtained in |14| . 
So it is not strange that the (|28|l should be identical to another nonspectral 
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equation. For completeness, the argument leading to that nonspectral equation 
is outlined. One notices that (|28l) is quadratic in terms of zi, and becomes 
identity when zi — Z2 or zi — z^. So (|28|l can be written as 

(zi - Z2){zi - Z3)Q{Z2, Z3) = 0, (31) 

which means (|28|l is equivalent to Q = 0. To obtain Q, one can simply put zi — 
in H28|) . Doing this, and inverting both sides, another equation is obtained which 
is quadratic in terms of Z3. Again it is seen that for Z3 = and Z3 — Z2, this 
equation becomes identity. So one can write this equation as 

23(^3- 22)0(22) = 0, (32) 

which is equivalent to Q = 0. To find Q, one simply writes the coefficient of 
z| in the inverted equation. One arrives at an equation containing only Z2. 
This in turn, can be converted to an expression quadratic in terms of 22- The 
coefficients of and z| of this equation are identities, while the coefficient of 
Z2 gives 

^2,3 bi,2ib2,3 + bi,2) = (62,3 + 61,2)62,3 61,2, (33) 

or 

62,3 [62,3, 61,2] = [62,3, 6l,2]6l,2, (34) 

which are the same as eqs. (47) and (48) in "T^, with 61,2 ^ 62,3, as expected. 



4 The conditional probability 



Assuming that the solvability condition H28|) . or equivalently H34|l . is satisfied. 
One can determine the conditional probability (or the propagator): 



U{x;t\y-0) 



i2n) 



N 



-wv 



,ia{p)-x tE(p) 



(35) 



where the integration region for each pi is [0, 2tt], and Ai = 1. The singularity 
in Aa- is removed by setting pj pj — ie, where the limit e ^ 0+ is meant. 
Note that as the elements of b are nonnegative and b satisfies (|10|l , the absolute 
value of the eigenvalues of b is not greater than 1. 

For the two-particle sector, there is only one matrix (b) in the expression of 
U. So, it can be treated as a c number. Using a calculation similar to what has 
been done in jT^ and one arrives at 



I7(f ; t\y; 0) = e" 



[xi - yi)! ix2 - y2)! 



y 

^ {.X2 - yi)\ (x 



y2-i)\ 



{x2 - vi)b 



t 



(36) 



To investigate the large-time behavior of the propagator, it is useful to de- 
compose the vector space on which b acts, into two subspaces invariant under 
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the action of b: the first subspace corresponding to eigenvalues with modulus 
one, the second corresponding to eigenvalues with modulus less than one. This 
is done introducing two projections Q and R, satisfying 



Q + R=l, 
QR = RQ = 0, 
[b, Q]-[6,i?] = 0. 



(37) 



Q projects on the first subspace, and R projects on the second. One can now 
multiply U hy 1 — Q + R. In the term multiplied by R, one can treat 6 as a 
number with modulus less than 1. But if the modulus of b is less than 1, then 
the integrand in H35|l is nonsingular and it is readily seen that for large t, the 
leading term in the integral is independent of 6, and in fact equal to the value 
obtained with 6 = 0. So 



U{x;t\y;0) 



-2t 



[xi - yi)! - y2)! 



fX2-yi 



-2t 



1=0 



{x2 - yiY- {xi - 2/2 - IV- 

lX2-y2 



-1 



{x2 -Vi)b 



Q 



[xi - yi)! (x2 - y2)! 



E 



fC2-yi 



fxi-y2-l 



{x2 - yi)! (xi - y2 - 0^ 



-1 



{x2 -Vi)b 



R, 
(38) 



and for large times, 

the second term = J_|e-[(^i-?'i-0"+(:r2-y2-t)"]/(2t) 
27rr 

_„-[{xi-y2-tf+(x2-yi-tf]/{2t) 



}R, t 



(39) 



So at large times, the second term tends to zero faster than t^^ , and the leading 
term in the conditional probability, which is of the order t^^ , does not involve 
the second term. 

If the only eigenvalue of b with modulus 1 is 1, then U has a simple behavior 
for t — > oo. In this case, bQ = Q, and one can simplify U to find 



[/(f;t|y;0) 



^xi-yi 



fX2~y2 



{xi 

OO 

E 

t - 



' yiV- ix2 - y2)! 

1X2 — yi j-xi—y2 — l 



{x2 - yi)! {Xl 

OO. 



y2-i)\ 



-1 



X2 - Vl 
t 



Q, 



(40) 



This is simply the propagator corresponding to a single-species drop-push sys- 
tem, multiplied by Q. 
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5 Dynamics of the particle-type 



For simplicity, let's continue with the two-particle sector. The probability that 
the first particle be of type Aa and the second particle be of type A^, regardless 
of their positions, is 

P«/5(<)= ^>/3(:z:i,X2;0, (41) 

where the primed summation means that the summation is on the physical 
region (xi < X2). Differentiating this, one arrives at 

P"/3(t) = ^ - 1, xr, t) + P"^(a;i, X2 -l;t)- 2P"'5(xi, t)], 

Xi ,X2 

= E(^"^- - - 1, a;; t), (42) 

X 

where in the last inequality, the boundary condition @ has been used. It is 
seen that the evolution equation of the particle-type is not closed; it involves 
the probability of finding different types of the particles in adjacent sites. 

However, the complete conditional probability for large times has the simple 
form (|40(l . For that form, the summation in H41|l is readily done, and one arrives 
at ^ 

J2 C/"''M-(2';t|y;0) = Q"V' t ^ ^- (43) 

a: — 1,2:2 

Here the fact has been used that the multiplier of Q in |@J is simply the propa- 
gator of the single-species drop-push model, and its summation on the physical 
region results in 1. From this, it is seen that the large-time probability of 
particle- types, depends only on the initial types of the particles, and not on 
their initial positions. This is of course true, when the only eigenvalues of b 
with modulus 1 is 1, the condition for 1)401) to hold. If moreover, this eigen- 
value is nondegenerate, then the large-time probability of particle-types is even 
independent of the initial particle types. In this case, Q would be written like 

Q"V = 9"^5^5,, (44) 

from which 

lim P"'3(t) = q"'^. (45) 

Here, q is the eigenvector of b with eigenvalue 1, normahzed as 

Sas^g"''-!. (46) 

It is seen that in this case, the large time probability of the particle types 
depends only on their interaction when they are adjacent. 
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6 Concluding remarks 



It was seen that a special class of multi-species drop-push models are solvable in 
the sense of the Bethe-ansatz. The condition corresponding to this solvability, 
resembles very much to what obtained in for the solvability of multi-species 
asymmetric exclusion processes. This is not accidental, since the behaviors of 
the drop-push model and the asymmetric exclusion model on continuum are 
related to each other: using a Galilean transformation, one can transform a 
drop- push model in which particles diffuse to the right, to an exclusion model 
in which particles diffuse to the left, TT, 12 . So the results obtained in |14) . 
regarding the solvability, with minor modifications can be used here. 
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